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A THEORETICAL ANALYSIS OF AIRPLANE LONGITUDINAL STABILITY 


AND CONTROL AS AFFECTED BY WIND SHEAR 

Windsor L. Sherman 
Langley Research Center 


SUMMARY 

Wind shear, the variation of horizontal wind velocity with altitude, has 
been a causative factor in several airplane accidents and may have been a con- 
tributing factor in other accidents. The longitudinal equations of motion with 
wind shear terms were used to analyze the stability and motions of a jet trans- 
port. A positive wind shear gives a decreasing head wind or changes a head wind 
into a tail wind. A negative wind shear gives a decreasing tail wind or changes 
a tail wind into a head wind. It was found that wind shear had very little 
effect on the short period mode and that negative wind shear, although it 
affected the phugoid, did not cause stability problems. On the other hand, it 
was found that positive wind shear can cause the phugoid to become aperiodic and 
unstable. In this case, a stability boundary for the phugoid was found that is 
valid for most aircraft at all flight speeds. Calculations of aircraft motions 
confirmed the results of the stability analysis. 

It was found that a flight-path control automatic pilot and an airspeed 
control system provide good control in all types of wind shear. Appendixes give 
equations of motion that include the effects of downdrafts and updrafts and 
extend the longitudinal equations of motion for shear to six degrees of freedom. 


INTRODUCTION 

Wind shear, that is, the change in the horizontal components of the wind 
with altitude, has been a causative agent in several airplane crashes that 
occurred during final approach*. (See refs. 1 and 2.) In addition to these 
accidents which caused a loss of life of 2^6 persons and injuries to many more 
people, wind shear may be an unidentified factor in many more accidents that 
occur during final approach and landing. Two of these accidents, the crash of 
Iberia Airlines, Flight 933, on December 17, 1973, at Logan International Airport 
(ref. 2), and the crash of Eastern Airlines, Flight 66, at John F. Kennedy Inter- 
national Airport on June 24, 1975 (ref. 3), are of particular interest because 
the wind shear profiles for these accidents are available; thus, an analytical 
investigation of airplane stability and control in accident-causing wind shears 
is possible. 

Considerable work has been done on the effects of wind shear (refs. 4 
to 7 ) , and it has concentrated on the effect on touchdown conditions . Gera 
(ref. 7) has formulated the longitudinal equations of motion of the airplane 
so that the dynamic effects of wind shear are considered . In reference 7 , the 


major concern was the effect on landing conditions and some stability considera- 
tions were included. 

The present paper adopts Gera's formulation of the problem as a starting 
point for the investigation of the longitudinal stability of the airplane, and 
the use of automatic control systems to improve the longitudinal characteris- 
tic3. The results presented are not as restrictive as would generally be sup- 
posed because the ainalysis of the accidents reported in references 2 and 3 indi- 
cated little or no effect on the lateral modes of motion. In addition to the 
work on longitudinal stability, a formulation of the six-degree-of-freedom equa- 
tions of motion based on Gera's approach is given. 

The airplane equations of motion used in the analysis are presented in 
appendix A. The airplane characteristics and flight condition are given in 
appendix B. Phugoid stability in wind shear is discussed in appendix C. 


SYMBOLS 

The SI system of units is used throughout this paper. All angles are in 
radians. 


c 

mean aerodynamic chord. 

m 

D 

distance along approach path, m; or d/dt 

E(D) 

numerator of airplane transfer function 

-> 

F 

force vector 


Fl »F2,F3 

forces along xi , X2, 

and X3, N 


forces along y-| 

and y3 axes, N 

g 

a 

acceleration due to gravity, m-S€C" 

H(D) 

denominator of airplane 

transfer function 

h 

altitude, m 


II 

moment of inertia about 

X-) , kg-m^ 

l2 

moment of inertia about 

X2 and y2, kg-m 

I3 

moment of inertia about 

, '3- 

X3, kg-m- 


i,j indices, 1,2, 3 

k,k-|,...k 7 gains 
9 , characteristic length, m 


2 



direction cosines 

= 1_ 35s, "'-'-sec-' 

I2 30 


I2 

9u 

1 

3M2 

I2 

da 

1 

3M2 

I2 

da 

1 

3M2 

I2 

36e 


M-|,M2,M3 moments about X'], X2, and X3, N-m 


m 

P 

P»q»r 

S 

s 

Ti/2 


mass , kg 
period, sec 

angular velocities about x-j , X2, and X3, rad-sec”^ 

wing area, m^ 

Laplace operator 

time to damp to half amplitude, sec 


Tdouble time to double amplitude, sec 
t time , sec 

Uq steady-state part of in stability axes, m-sec“^ 

u changing part of in stability axes, m-sec"' 

V airplane speed along stability axes, m-sec“^ 

airplane velocity- vector, m-sec”' 

Vg vector acceleration in Earth-fixed axes, m-sec“^ 

resultant velocity vector, m-sec~^ 
wind velocity vector, m-sec“^ 

^G, 1 »^G, 2 »^G ,3 components of in Earth- fixed axes, m-sec“^ 

vi,V2,V3 components of along x-j , X2, and X3, m-sec“^ 
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^w, 1 »^w, 2 >^w ,3 components of along X-| , X2, and X3, _rn -sec~l 

^w,1»^^w,2 wind gradients with altitude, sec“^ 

wind gradient along flight path, sec“^ 
initial wind speed, m-seol 

F„ 1 ^_i 




'■a 




= 1 -ZjJ., sec" 
m du 

= J. ^y» ^ , m-rad“^-sec“2 
m da 

= 1 ^y,> ^ , m-rad“^-sec“^ 
m 8 60 


X-|,X2,X3 Earth-fixed axes 

xj,X2,X3 nonrotating axes that translate with air and are always parallel to 
X-| , X2, and X3 

x-|,X2,X3 airplane principal body axes 
yi>y2>y3 stability axes 


Zq 

= m-3^’ 

m-rad“^-sec“^ 

Zu 

= 1 ^y,3, 
m 9u 

sec“^ 

Za 

= 1 Im, 

m da 

m-rad“^-sec“2 

z<i 

= 1 ^y,3, 
m dd ’ 

m-rad“^-sec“^ 

2«e 

= l!l3U3, 

m d6g 

m-rad“^-sec“2 


a 

«0 

Aa 

e 

r 


angle of attack, Oq + Aa, rad 
steady-state angle of attack, rad 
change in angle of attack, rad 
angle of sideslip, rad 
flight-path angle, Fq + Y> nad 
flight-path command - 
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ro steady-state flight-path angle, rad 

Y change in flight-path angle, rad 

6@ elevator deflection, rad 

flap deflection, rad 
Cp phugoid damping ratio 

0 Euler pitch angle, rad 

p air density, kg-m“3 

ax total shear parameter 

Oy vertical wind shear parameter 

a„ downdraft-updraft parameter 

T servo time constant in flight-path control system - 

Tg servo time constant in speed control system - ^ C 

Euler roll angle, rad 
Euler yaw angle, rad 

w airplane angular velocity vector, rad-sec”^ 

0)p .phugoid natural frequency, rad - Sec" * 

Subscript: 

1 initial 

Dots over symbols indicate differentiation with respect to time. An arrow 
over a symbol denotes a vector. 


ANALYSIS 

Work in wind shear is somewhat complicated by the lack of positive defini- 
tions; therefore, it is necessary to define exactly what is meant by the termi- 
nology used herein. Vertical wind shear is the change of the horizontal compo- 
nents of the wind with altitude and is referred to as wind shear. Head winds are 
negative and tail winds are positive. A positive wind shear gives a decreasing 
head wind or changes a head wind into a tall wind. A negative wind shear gives a 
decreasing tail wind or changes a tail wind into a head wind. This definition is 
maintained whether the airplane is climbing or diving. 
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The variation of wind speed with altitude that occurred on the approach 
path at John F. Kennedy International Airport at New York on June 2M, 1975, is 
shown in figure 1 . The data given in figure 1 were recorded about the time of 
the crash of Eastern Airlines Flight Number 66 and, because the wind gradient 
v^ 1 exceeds 0.1 sec“^, these data represent a severe wind shear. In figure 1, 
head winds are negative and tail winds are positive. From 178 m to 116 m there 
are rapidly changing head winds. At an altitude of 116 m, a well-defined shear 
starts and continues down to 50 m. Table I shows the Vw,1 that existed over 
various altitude intervals. (See fig. 1.) 

TABLE I.- WIND SHEAR GRADIENTS 
[Taken from fig. lj 


Altitude interval, m 

Vw,1. sec-1 

116 to 102 

0.10 

102 to 80 

.35 

80 to 68 

.09 

68 to 50 

.27 

Average over altitude range . . . 

0.20 


A useful parameter in the analysis of wind shear is the wind shear parame- 
ter Oy. This parameter which comes from the linearization of the equations of 
motion (see appendix A) is defined as 

(O 

g 

When equation (1) is applied to the gradient given in table I, it is found that 
for a typical commercial jet with an approach speed of 77.12 m-sec“^ , 
varies from 0.79 to 2.75. 

When a vertical wind shear, that is, the variation of horizontal wind speed 
with altitude, occurs in a thunderstorm, it is likely to be accompanied by 
updrafts and downdrafts. Severe downdrafts are sometimes referred to as down- 
blasts. The updrafts and downdrafts that occurred along the approach path at 
Kennedy Airport at the same time as the vertical shear shown in figure 1 are pre- 
sented in figure 2. The corresponding parameter for an updraft or downdraft is 

a„ = (2) 

g 

where ^W, 3 is the gradient. Figure 2 was used to determine gradients for the 
downdrafts and updrafts and the gradients were found to vary between 0.05 sec“^ 
to -0.06 sec"^ . The corresponding values of were 0.31 and -0.47. The data 

obtained from figures 1 and 2 and reference 8 were used in the analysis of the 
longitudinal stability of a large four-engine jet transport. The characteris- 
tics of this airplane and the flight conditions used are given in appendix B. 
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The changing wind speed imparts an acceleration to the reference system 
(that is, the nonrotating axes that translate with the airplane); thus, this 
axis system is.no longer inertial and Earth-fixed axes must be used. The total 
acceleration Vg is given by 

Ve = Va + (3) 

A set of longitudinal equations of motion based on equation (3) are given 
in reference 7. These equations take into account vertical wind shear, but 
downdrafts and updrafts are not considered. In appendix A, the equations pre- 
sented in reference 7 are extended to include the effects of downdrafts and 
updrafts. Equation (A18) was used to study the longitudinal stability of a 
large jet transport. The total shear parameter is • introduced in equa- 

tion (A18). This parameter is defined as 

Ot = % + a„ (4) 

In order to obtain a clear understanding of the effects of vertical wind shear, 
the analysis was made with = 0, and the effects of Oy were considered 
after the vertical shear study. 

Equation (A18) is very interesting as the presence of Ox and in the 

dynamic terms indicates an interaction between the airplane and its environment. 
This equation was used to make a stability analysis of the large four-engine jet 
transport described in appendix B. The analysis was made for flight-path angles 
for 0 to ±0.1745 rad in steps of 0.01745 rad. The range of values for Ox 
(-3.5 g Ox ^ 3.5) was based on the wind shear data presented in figure 1i 

The left-hand side of equation (A18) was used to obtain the characteristic 
equation of the airplane, a quartic, which was solved to determine the stabil- 
ity. The results obtained for Fq = -0.05236 rad are given in table II. Fig- 
ure 3, a root-locus plot based on table II, shows the motion of the roots as 
Ox changes. The data presented in table II and figure 3 show that negative 
wind has little effect on the airplane. In the case of positive shear, these 
data Indicated no appreciable effect on the phugoid until Ox exceeds 1.0. At 
Ox = 1.0, the phugoid has broken down into two aperiodic modes, one of which is 
unstable. The stable mode damps to half amplitude in about 54 sec, whereas the 
unstable mode takes 333 sec to double amplitude. As Ox increased above 1.0, 
the time to double amplitude decreased until it was 3*1 sec for Ox = 3.5. The 
results for the other flight-path angles used in the study were consistent with 
those obtained for Fq = -0.05236 rad. 

The breakdown of the phugoid oscillation into two aperiodic modes for 
Ox 1.0 was traced to the term g(cos Fq - Ox cos 2Fq) which changes sign at 
this value of Ox. The change in sign of this term is independent of the sign 
of Fq so the same type of instability will occur whether Ox ^ 1.0 for climb- 
ing and level flight. A more complete discussion of this term is given in appen 
dix C. The effect of wind shear on the short-period longitudinal modes was 
negligible. The short-period roots were periodic -0.70058998 ± 0.809484271 for 
Fq = du = 0 and -0.68874419 ± 0.78l30806i for Fq = -0.05236 and dy =,3-5. 
Increasing the flap deflection to 0.87266 rad produced very little or no change 
in these results. 
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TABLE II.- EFFECT OF POSITIVE AND NEGATIVE SHEAR 


ON PHUGOID MODE - BASIC AIRPLANE 
j^iSf = 0.4363 rad; a„ = O.Oj 


To, 

c^u 

Roots 

Ti/2» 

"^double > 

P, 

Wp 9 


rad 




sec 

sec 

sec 

rad 


0.0 

0.0 

- 0.002954 

±0. 140281 

234.59 


44.79 

0.14031 

0.021 

-.05236 

.0 

-.0052453 

±.140501 

132.09 


44.72 

.1406 

.037 



-.5 

-.0052994 

±.171471 

130.77 


36.64 

.1716 

.031 



-1 .0 

-.0054139 

±.197251 

128.00 


31.85 

.1973 

.027 



-1.5 

-.0055879 

±.219691 

124.02 


28.60 

.2198 

.025 



-2.0 

-.0058200 

±.239741 

119.07 


26.21 

.2398 

.024 



-2.5 

-.0061076 

±.257971 

113.47 


24.36 

.2580 

.024 



-3.0 

-.0064496 

±.274751 

107.45 


22.87 

.2748 

.023 



-3.5 

-.0068442 

±.290321 

101.25 


21 .64 

.2904 

.023 



.5 

-.0052567 

±.0996191 

131 .83 


63.07 

.9962 




1 .0 

-.012747 

.0020821 


332.84 






1.5 

-.10647 

.095524 


7.25 






2.0 

-.14893 

.13756 


5.04 






2.5 

-.18207 

.17013 


4.07 






3.0 

-.21051 

.19785 


3.50 






3.5 

-.23600 

.22249 

■ 

3.11 





Phugoid Stability Boundary 

As pointed out in the previous section, a study of the calculations 
that produced table II showed that unstable conditions occur when the term 
g(cos Fq - O’j; cos 2Fq) in equation (AI 8 ) changes sign. The only possible way 
in which this term can change sign is for Oj to become greater than a critical 
value. For Fq o 0, the critical value of (Jf is 1.0. Additional calculations 
showed that for the normal range of flight-path angles, the difference between 
the critical value of (J-p and 1.0 was negligible. These results lead to a sta- 
bility boundary for the phugoid in wind shears that is independent of airplane 
characteristics. (See appendix C.) This boundary is obtained by setting equa- 
tion (1 ) equal to 1.0 and solving for as a function of the speed. This 

boundary is shown in figure 4(a) . The approach speed bands for general aviation 
and jet transports have been marked in figure 4(a). A glance at figure 4(a) 
reveals that the higher the speed of the airplane, the smaller the wind gradient 
for the onset of unstable conditions. Because of the slower approach speeds, 
general aviation aircraft can tolerate a higher gradient, before becoming unsta- 
ble, than transport aircraft. As there were no airplane characteristics used to 
determine the curve shown in figure 4(a) , it is necessary to assume a specific 
airplane in order to determine stability conditions around the stability 
boundary. 
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A segment of the curve from 50 m-seo"^ to 100 m-sec"^ is shown in fig- 
ure 4(b). Times to double amplitude (above the curve) and times to half ampli- 
tude (below the curve) have been calculated for the airplane described in appen- 
dix B at three specific speeds, 67 m-sec“^ , 77 m-sec"^ , and 87 m-sec“^ , as a 
function of wind gradient. This speed range covers the approach speeds for jet 
transports. As would be expected, the time to double amplitude decreases as the 
gradient increases. As the speed is increased for a given value of v^^l that 
is near the stability boundary curve, changes in stability are most drastic. 

For ^w,1 = 0 .15 the time to double the amplitude decreases from 35.1 sec at 
67 m-sec“^ to 8.9 sec at 87 m-sec“^. At v^ 1 = 0.125, the phugoid changes frofli 
stable to unstable over this speed range. T^ese results are interesting as they 
indicate that the phugoid stability can be helped by slowing the airplane down, 
a practice not recommended because the airplane is approaching its stall speed. 
On the other hand, increases in speed, particularly for gradients close to those 
for aj = 1*0, can cause an unstable condition to develop. 


Use of Automatic Pilots 

Automatic pilots offer a means of improving the airplane characteristics 
when it is flying through wind shear. There are two basic types of autopilot 
that can be used to control longitudinal motion. These are a 0 autopilot that 
controls pitch attitude and a F autopilot that controls the flight path of 
the airplane. The 6 autopilot was not considered for the following reasons: 

(1) During landing approach, a constant flight path is desired and a 6 
autopilot provides only indirect control of the flight path. 

(2) A 0 autopilot places a constraint on the pitch attitude of the air- 
plane and forces all flight-path variations into the angle of attack. The vari- 
ations in a can have an adverse effect on lift and may even cause the airplane 
to stall. 

A block diagram of the F autopilot used in this study is presented in 
figure 5. 

There are four gains and a servo time constant to be assigned values. The 
values used are 

k = 0.8 

ki = 5.0 

k2 = 1.0 

k3 = 7.7 

T = 0.0357 

and were taken from the successful F autopilot used in reference 9. By using 
these gains and the data presented in appendix B, the roots of the characteris- 
tic equations were determined for the same range of Fq and Oj as the basic 
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airplane or in table II. The characteristic equation is a six-degree equation. 
Again, these calculations were made with = 0.0. The negative wind shear 
case will be considered first. A typical set of roots when Fg = -0.05236 rad 
and 0T = -1.00 are . 

0.012992936 

-0.27539026 ± 0. 145479611 
-0.46845553 ± 2.1273381 
-27.956636 

The two periodic roots damp to half amplitude in 2.52 and 1.48 sec, respec- 
tively. The stable aperiodic root damps to half amplitude in about 0.025 sec 
and it takes the unstable aperiodic root about 53 sec to double amplitude. This 
root pattern was constaint throughout the negative shear calculations. The unsta- 
ble condition comes from the term 

^ sin 2Fo - Xuj (-Za) - ^ sin2 Foj(-Xa) 

which is the constant term in the characteristic equation of the airplane auto- 
pilot system with and Zg^ set equal to zero. The foregoing results were 

consistent over the range of negative shear and flight path considered in the 
study . 

When the wind shear was positive, the same root pattern was found for small 
values of Fg and ^0 = -0.05235 and Of □ 1.0, the roots are 

0.00099216978 

-0.27327956 ± 0. 120537261 

-0.45804227 ± 2.1l4l669i 

-27.956737 

which, with the exception of the unstable root, are very close to the roots 
obtained for negative shear. When the shear factor was increased to 1.5, the 
airplane autopilot combination became stable and the roots are 

-0.0024078146 

-0.27256525 ± 0.112939991 

-0.455542571 ± 2.11085721 

-27.956762 

This small change in stability occurred because the term (Cxg/2Ug) sin 2Fg 
has become large enough to change when combined with -X^, the sign of the con- 




(5) 
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stant term (expression (5)); thus, a stable system is obtained. The foregoing 
calculations were made with Xg and Zg equal to zero. Both cases were 
rerun with estimated values for®these parameters and no significant changes were 
found in the results. The results obtained for positive shear were most inter- 
esting as a wind shear that destabilized the phugoid for the airplane alone now 
exerted a stabilizing effect for the airplane autopilot combination. 

Because the unstable mode in the airplane-autopilot system occurred for 
Fq = ax = 0.0, this condition indicated that the source of the problem was the 
airplane-autopilot system and not the wind shear. 

For Xg = Zg = 0 and ax = Tq = 0, the constant term in the characteristic 
equation of the airplane autopilot system is (X^Z^j - ZuXo[)Mg kk2 and the unsta- 
ble condition occurs for the configuration used in this study, because Z^X(^ is 
greater than X^Z^j. An increase in XyZ^j will restore stability. For the air- 
plane used, this increase in Xu?a can be obtained by increasing the flap 
deflection from 0.4363 rad to 0.87266 rad. This doubling of the flap deflection 
gives XyZ^j = 2.40642 and ZyX^ = 1.88339 and results in a positive differ- 
ence and a completely stable airplane autopilot combination. The roots of the 
characteristic equation for Fq = 0.05235 rad and a negative shear, factor 
aj = -1.0 are 

-0.0028802057 

-0.27307666 ± 0.15105564i 
-0.48581279 ± 2.1097746i 

-27.949297 

The instability just discussed arises because of the flight-path constraint 
imposed on the airplane by the flight-path control system. This type of insta- 
bility, stability under partial constraint, is discussed in detail in refer- 
ences 10 and 11. The method used to eliminate it, that is, increasing flap 
deflection, is the same as that suggested in reference 12. 

V/hen the shear factor increased to -1.5, an unstable root again appeared in 
the characteristic equation. The shortest time to double amplitude was approxi- 
mately 124 sec, a length of time that should pose no problems of control. 

Depending on the value of Fq, the transition from stable to unstable conditions 
occurs at different values of aj; the smaller Fq, the larger the ax required 
to cause an instability. When the shear factor was positive, ax = 1.0; for the 
same Fq the roots of the characteristic equation are 

- 0.016587880 

-0.26976410 ± 0.12312273i 
-0.47574423 ± 2.0964l92i 

-27.949405 
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Over the range of Fq and Oj considered in the investigation, the airplane 
autopilot system remained stable; however, as both Fq and increased, the 

size of the first real root increased and decreased the time to damp to half 
amplitude. 


Updrafts and Downdrafts 

*• 

As pointed out at the beginning of the analysis section, updrafts and down- 
drafts can occur (see fig. 2) in conjunction with vertical wind shear. Like 
vertical wind shear, gradients can occur in upward and downward vertical veloci- 
ties. An analysis of the measured data presented in figure 2 showed gradients 
from 0.06 sec"^ to 0.05 sec“^. For these gradients, the downdraft parameter 
lies in the range -0.47 ^ ^ 0.31. 

When Oy is included in the equations of motion, it appears as an additive 
factor to (7u> *^T = “^u + '^w» except in the coefficient of the u term of the 
second equation of equation (Al8). In this coefficient the term (g/Uo)a^^ 
appears with no corresponding term. The effect of this term on the results 

is negligible. For instance, let = 2.0 and a„ = 0.4. Calculations showed 

that with the (g/Uo)d^^ omitted, = 2.4 and the predicted time to double 

amplitude was 4.26 sec. When (g/Ug) ^w was included, the predicted time to 
double amplitude was 3-82 sec, a change of 0.44 sec which is negligible. Calcu- 
lations using the range of du and noted in this report should show simi- 
lar results . It appears that unless becomes very large , an excellent 

approximation of the stability can be obtained if the term (g/Uo)0;^ is not 
considered . 


The only change, outside of the summation of and a„, in the quadratic 

approximation for the phugoid mode is in the expression for 03 p 2 which becomes 



Xy(sin Fq - dj sin 2Fg) - 


(Zu - L 


(cos Fq - d-j COS 2Fq) 


( 6 ) 


where d>p = d^ + d„. If (g/UQ)d„ > Zy the stability is still determined by 
the sign of cos Fq - dx cos 2 Fq. If, however, (g/Uo)d„ < Zy, sign changes 
will occur that will affect the prediction of stability. Since (g/Uo)d„ = 2» 

this means that Vw, 3 would have to be less than -0.3. For the severe wind con- 
dition used in this report v ^^3 did not exceed ±0.1 sec”^. As Zu varies from 
-0.15 to -0.40 for most airplanes, it is concluded that only in rare Instances 
will the term (g/Uo)d„ influence the stability calculated from equation (A18) 
and the phugoid approximation. 


The phugoid stability boundary (fig. 4(a)) is based on the zeros of 
in that Up^ will change sign for dx > 1*0. This boundary was recalculated by 
using equation (6) and it was found that over the range of d^, values consid- 
ered in this paper (which are for severe weather conditions), (g/Uo)d^ did not 
change the result that over a normal range of flight-path angles, the sign of 
0)p2 is dependent on dx. Thus, for most the phugoid stability boundary of 

figure 4(a) is valid. 
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DISCUSSION 


In the preceding section, an analysis of the longitudinal stability of a 
large transport airplane in wind shear was made. The results of this analysis 
definitely show an interaction of the airplane and its environment. The inter- 
action comes from accelerations produced by the wind gradient. The gradient is 
the important wind factor, not the wind speed. These facts made the results 
obtained difficult to interpret because the classical stability conditions are 
concerned with the stability about a steady-state condition and wind shear is a 
transient phenomenon. The equations used for the analysis give the stability 
over regions of constant wind shear. 

The results of the stability analysis show that shear, particularly posi- 
tive shear, can have a drastic effect on the stability of an airplane. In order 
for the pilot to take action, it is necessary for him to know whether a wind 
shear is present and what type (that is, positive or negative) wind shear is 
affecting the airplane. References 13 and 14 attack the problem of measuring 
the total energy rate of the airplane and describe devices for its measurement , 
the one described in reference 13 being particularly simple and effective and 
suitable for general aviation aircraft. Reference l4 goes one step beyond ref- 
erence 13 and shows how the change in total energy can be used to detect wind 
shears . 


Control-Fixed Aircraft Mo'tions in Wind Shear 

In order to obtain better insight into the meaning of the results obtained 
from the analysis made with the stability equations, the six-degree-of-freedom 
equations (see appendix A) were reduced to the three-degree-of -freedom equations 
of longitudinal motion with = 0.0 and programed for digital computations of 
aircraft motions. 

In all cases, including those with head winds and tail winds, the airplane 
was trained on a -0.05236 rad flight path at an altitude of 130 m. Thus, the 
only disturbances that occur during the flight are the vertical wind shears that 
start at an altitude of 106 m and end at an altitude of 50 m. 

Plots of altitude, pitch ahgle, and flight path against D are shown in 
figure 6. The distance along the approach path* D was used Instead of time for 
the plotting since this distance relates the motions to the space position of 
the airplane. All the motions shown in figure 6 are for fixed controls. 

Curve A, which is for zero initial wind speed, = 0, and Oq; = 0.0, is taken 
as the reference case. A 6.1-m-sec“^ head wind (curve B) causes the airplane to 
undershoot and the glide slope becomes a little steeper; the pitch angle is 
larger, but remains relatively constant. 

When a tail wind of 6.1 m-sec"^ is encountered (curve C) , the airplane 
floats and the flight path is not as steep as the reference case.. As would be 
expected, the magnitude of the pitch angle is less than was found for the refer- 
ence condition of curve A. In all three cases the curves are smooth and do not 
exhibit oscillatory characteristics. The remaining cases (curves D to G) show 
the effects of shear when combined with a head wind or tail wind. Curves D and 
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E are for head winds and curves F and G are for tail winds. In all cases the 
shear started at an altitude of 106 m and continued until the airplane reached 
an altitude of 50 m where the airplane flew out of the shear. In the case of 
strong positive shear, Ox w 2.0, the altitude, 0, and y motions are charac- 
teristic of a system with aperiodic roots, one of which is unstable, and the 
airplane Impacts about 1500 m short of the undisturbed impact' point. From the 
onset of the shear, only 7 sec elapsed before impact. When the shear was 
decreased, Oj = 0.5, the characteristic of the traces changed, and the airplane 
was just starting to level off in what appears to be the start of an oscillating 
response when impact occurred about 800 m short of the undisturbed touchdown 
point. The stability analysis for this case predicted an oscillating response 
and the impact condition is consistent with long time to damp to half amplitude, 
about 65 sec, predicted by the stability analysis. For curves F and G, W^ was 
changed to 6.10 m-sec"^ and the shear parameters aj were -0.5 and -2.0. As 
would be expected, the airplane floated and a large overshoot occurred. As pre- 
dicted by the stability, the response is oscillating but damped, the oscilla- 
tions being induced by the wind shear which excited the phugoid mode. As a 
result of these motion studies, the results that were predictable from the sta- 
bility analysis, that positive shears are most dangerous to aircraft and that 
negative shears are much less so, have been confirmed. 

Many airplanes in the past have been statically unstable in the power 
approach condition. When such an airplane is allowed to diverge downward with 
fixed controls, the airspeed increases and the angle of attack decreases. In 
the case of a down divergence in positive wind shear, however, the airspeed 
decreases and the angle of attack increases, and the normal recovery technique 
of the use of up elevator to recover results in stall and further loss of alti- 
tude. The divergence due to wind shear, although it affects the constant term of 
the stability equation the same as static instability, would result in entirely 
different piloting problems. 

Motion results were also obtained for a flap setting at 0.87266 rad instead 
of 0.43633 rad with little change in the responses of the airplane. The effects 
of downdrafts and updrafts were found to affect the airplane in the same manner 
as positive and negative shears, respectively; thus, the shear parameter need 
only be adjusted to account for this effect. 

The usual method for simulating wind, turbulence, and other wind effects., 
such as shear, in airplane simulations is to include the effects in the resul- 
tant velocity, vector used in calculating the airplane forces and moments. Terms 
such as the terms in equations (A4) to (A6) are generally not included in 

the equations of motion. These terms that appear in equations (A4) to (A6) 

or their equivalent were included in all calculations made in this study. Fig- 
ures 7 and 8 show the effect of neglecting these terms for positive and pegatlve 
shear. Although the general character of the motions is the same, the 
terms have a significant effect on the airplane motions. 


Use of Autopilots To Control Airplane in Wind Shear 

The results of the stability analysis indicated that an automatic flight- 
path control system provided an increase in the stability of the airplane. The 
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motion study showed that changes in airspeed as large as 23 m-sec”^ were occur- 
ring for the shear passage, which implies that an airspeed control is necessary 
so that the lift of the airplane is maintained as constant as possible. The 
flight-path and speed-control systems were checked separately before being com- 
bined. The flight-path control system (fig. 5) was added to the simulation and 
the results showed large changes in angle of attack that caused the airplane to 
stall. Changing the f feedback to a q feedback provided increased stability 
and eliminated the large excursions in angle of attack. All results presented 
in the section are for the flight-path control system with a q pitch rate feed- 
back. The speed-control system was set up to keep the resultant speed, that is, 
the sum of wind-speed vector and ground-speed vector constant. The resultant 
speed (or airspeed) is available from the air data computer. As shown in fig- 
ure 9, the engine was simulated by a first-order servo. The time constant of 
this servo was selected so that 90 percent of the command thrust was achieved in 
5.6 sec and a xg of 2.5 resulted. This time delay matches well with the 
engine response data given in figure 28 of reference 15 for engine response to 
throttle. A Xg of 0.1 was also used to simulate the thrust modulators of ref- 
erence 15. The values of Xg given were assumed to simulate the extremes of 
engine response. The results of these calculations are shown in figures 10 and 
1 1 . The flight-path and airspeed control systems successfully guided the air- 
plane through a wind shear where Oj = ±2.0. The positive shear case shows more 
sensitivity to the engine time constant than the negative shear case. The speed 
control system, through its input command, attempts to hold airspeed constant. 

In combination with the acceleration feedback, this control system also attempts 
to hold acceleration to zero. This procedure is roughly the equivalent of 
attempting to hold jifj} • ^gl/g constant, which is part of the energy rate equa- 
tion given in reference 11. Good results were obtained by yse of this autopi- 
lot. It appears that a critical term to control is I^r * V^j/g and that more 
work should be done on automatic control systems that directly control this 
term. 


CONCLUDING REMARKS 

The effect of the phenomenon of wind shear, that is, the variation of 
horizontal winds with altitude, on the stability and motions of a large jet 
transport has been investigated. In the case of the phugoid mode negative wind 
shear, a shear that changes tail winds to head winds, the phugoid remained peri- 
odic and stable, although the time to damp to half amplitude decreased as the 
shear gradient increased. Positive wind shear, a shear that changes a head wind 
into a tail wind, caused the phugoid to become unstable for values of the shear 
parameter greater than 1.0. The motion was aperiodic and the time to double 
amplitude decreased as the shear parameter increases. In the case of positive 
wind shear, a stability boundary for the phugoid, valid for all aircraft, was 
found as a function of the wind gradient and approach speed. The stability 
boundary is independent of airplane characteristics and, apparently, an attempt 
to hold airspeed constant in wind shear conditions is the best action to avoid 
changes in airplane stability. 


The effect of wind on airplane motions confirmed the effects predictable 
from the stability analysis. More study is needed of autopilots for the ade- 
quate control during a wind shear encounter. The results Indicate that it is 
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important for the pilot to know whether a wind shear exists and whether it is 
positive or negative. 

Downdrafts and updrafts were considered as well as horizontal wind shears. 
It was found that downdrafts and updrafts have the same effect as wind shear and 
could be accounted for by modification of the shear parameter. 

A stability analysis of the airplane equipped with a flight-path control 
autopilot showed that the combination was a stable system. When the shear was 
negative, increasing the flap deflection helped system stability. Motions 
obtained with a flight-path control with a pitch-rate feedback and a speed- 
control system that held constant speed with respect to the air mass showed 
that with these control systems, an airplane could fly through vertical wind 
shears with no adverse effects. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton , VA 23665 
April 27 , 1977 
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APPENDIX A 


AIRPLANE EQUATIONS OF MOTION USED IN THIS STUDY 


The linear accelerations of an airplsine in a moving axis system are given 
by the vector equation 

Ve = ^(^A + + (3 X + ^w) 

at 

and the corresponding force equations are 

1_(Va + V„) + cS X = I (A1) 

dt m 


Equation (A1) is general and may be resolved into any desired coordinate 
system. Because the use of principal body axes simplifies the resulting equa- 
tion,, these axes were used for the six-degree-of-freedom equations. The wind 
vector is defined in Earth-fixed axes, X-| , X2, and X3, head winds are 

negative, and is defined in principal body axes, x^i , X2, and xj. The 

X-| and x-| axes are positive in the direction of flight and the X3 and X3 
axes are positive downward. (See fig. 12.) 


In order to obtain the force equations in the principal body axes, it is 
necessary to transform the wind from the fixed-axis system to the body-axis 
system. This is done by using the nine direction cosines obtained, from the 
standard 0, and 9 Euler transformatioft. The direction cosines Jlj[j are 
obtained from 




Olij)i 



(A2) 


where 



direction cosine 


rates. 


are given by 




J-il = r£-i2 - U^13 
ii2 = P^!'i3 - 




(A3) 


^13 = <1^11 - P%2j 

By using equations (A2) and (A3), equation (A1) may be written in component form 
as 


VI + Vwjilii + Vw,2^21 + + qv3 - rV2 = ^318 + ^ (A4) 

V2 + v„jJ,i2 + ^,2^22 + ^v,3h2 + - pvg = £328 + ^ (A5) 
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V3 + + v„, 2^23 + ^w,3^33 + P^2 - <1^1 


= A33S + 


h 

m 


(A6) 


The variables p, • q, and r are the angular velocities of the airplane in prin- 
cipal body axes and are obtained from the equations of angular ipotion which, in 
component form, are 


P + ^3 - ^2 qr = !!l 
II II 


Ii - lo Mp 

q + _! £ pr = _£ 

I 2 I 2 


(A7) 


(A8) 


r . ^2 - Ii pq = ^ 

I 3 I 3 


(A9) 


The aerodynamic forces and moments are functions of the resultant velocity 
given by 

= ^A - 

as are the parameters a, d, 3> and 3 . Because of the wind shear stability 
derivatives normally neglected (such as Cj)^, Cjj^, C 2 ^p, and C^qj) should 

be included in the formulation of the forces and moments . The wind vector Vw 
includes steady winds, turbulence, and wind shear. Head winds are negative. 


An examination of equations (A4) to (A6) indicated that in addition to the 
usual terms, the terms v^ 1 , v„j 2 » ^^w,3 required. These terms are 

wind acceleration terms an(i may arise from turbulence or wind shear. In models 
for turbulence, it is fairly easy to write the model so that the accelerations 
are available for equations (A4) to (A6). In the case of wind shears, particu- 
larly where recordings of actual wind shears are used, obtaining the required 
accelerations is not quite so straightforward. In order to have a consistent 
method of specifying wind shears, the following method is suggested. An acceler- 
ation may be written as a product of a velocity gradient taken over some charac- 
teristic length and the rate of change of this length. Thus, 


a - dv dil 
dil dt 


(A10) 


The length Jl over which the gradient is determined should be 30.48 m, a length 
that is consistent with the reporting of vertical wind shears. Thus, winds sepa- 
rated by 30.48 are subtracted and divided by 30.48 to obtain dv/dA which has 
the dimensions of sec“^ . Thus, wind accelerations for equations (A4), (A5), and 
(A6) are given by 


V 


w,1 


dvw,1 

dJl 


Vg,3 


(All) 


18 



APPENDIX A 


^w,2 = 


dv, 


w_j_2 _ 

dfi, 


(A12) 


'w,3 


dJl 


Vg ,1 


(A13) 


where dv„ i/dJl, dv„ 3 ^*^^ velocity gradients and usu- 
ally written as Vw, 1 » and v^^ 3 . 

Equations (A4), (A 6 ), and (A 8 ) are the longitudinal equations of motion and 
were used to calculate the airplane motions presented in this report. 


For the stability calculations presented in this paper, the longitudinal 
equations of motion were transformed to stability axes. (See fig. 12.) After 
making this transformation and substituting the trigonometric equivalents for 
the direction cosines and using the definition of F, F = 6 - a, the longitudi- 
nal equations of motion in stability axes are 


V - 


Vw, 1 sin r cos r - 3 cos 


F + g sin F = ^ 

m 


(A14) 


-VF - v^^iV sin^ F + v^^3V cos^ F 


g cos F = ^y>3 
m 


(A15) 


f + a = ^^2 (A16) 

I 2 

In these equations a single component vector in the y-] -direction 

in stability axes, and called V. Equations (A14) to (AI 6 ) were linearized by 
making the usual assumptions that 

V = Uq + u 


oi = ttQ + Aa 


F = Fo + Y 

where the zero subscripted terms are the steady state and u, Aa, and y are 
the perturbations from the steady state. For the linearized equations, it is 
convenient to use Oy and which are constants in place of v^^i and 

v^^ 3 . These parameters are defined as 

g 

= Pqv4,3 
8 
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The linearized equations of motion are 


sin 2 F 0 

dt l2Uo 2Uoj 


u + /gjcos Tq - (dy + cos 2roj\ Y = ^ 


sin^ Tq + cos^ Tqiu + 

Uo Uo 


(du + d„) sin 2 rq] ? Y = 


it 


> (A17) 


_^£_(Y + Aa) = ^ 
dt2 I 2 




The forces and moments were expanded and in matrix form these equations, 
writing a for Aa and o-p for Oy + cr^j become 


^ ^ ax Sin 2Fo - 


-Zy - i_(ax sin^ Tq - a^) 

Uo 




-Xa 


-(Z^ + “ 2q( 


^ - iHi * V|t - It. 


g(oos Tq - ax C03 2Tq) 


-(Oq + 2q)^ + g(sin Tg - ax sin 2rg) 


(jt2 dt 


*6, 






«e (A18) 


L J 


For 0 ^= 0 , these equations reduce to the equations given in reference 8 
and for d^ = d^ = 0, they become the usual longitudinal equations of motion. 
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AIRPLANE CHARACTERISTICS AND FLIGHT CONDITION 

The airplane used in this study is considered a typical narrow body modern 
jet transport airplane powered by four engines, each having approximately 
67 233 N of thrust. 

The dimensional and mass characteristics are 
c = 7.01 m 
S = 267.9 m2 
rr.= 90 909.1 kg 
l2 = 9 933 300 kg-m2 
p = 1.2929 kg-m-3 

The aerodynamic data for the stability axes, center of gravity at 0.25c, are 



0.43633 rad flap 

0.87266 rad flap 

Za 

-55.055 m-rad“^-sec”2 

-52.68 m-rad-^ -sec-2 

Zu 

-0.29024 sec- I 

- 0.29024 sec-’’ 

Zq 

-3.2708 m-rad-^-sec”* 

- 3.2708 m-rad-^ -sec-2 

Zd 

-1.0075 m-rad-^-sec”^ 

- 1.0075 m-rad-^-sec-^ 

2«e 

-2.63428 m-rad“^-sec-2 

- 2.63428 m-rad-^ -sec-2 

Xa 

- 5.9803 m-rad“^-sec-2 

-6.48907 m-rad-^ -sec-2 

Xu 

- 0.02385 sec-l 

-0.04568 sec-^ 

^«e 

-0 . 1 568 m-rad“ ^ -sec~2 

-0.1568 m-rad-^ -sec-2 

Ma 

- 0.809 rad“^-sec-2 

-0.8468 rad- ^ -sec-2 

Mq 

-0.51 3 rad” ^ -sec” ^ 

-0.5481 rad- ^ -sec- ^ 

% 

-0.175 rad-l-sec-1 

- 0.18778 rad- ■'-sec- "I 

Mu 

-0 . 00095 m-^-sec”^ 

-0.00095 m-^-sec-^ 

M6e 

-0.73733 rad- ■'-sec-2 

- 0.75038 rad- ■'-sec-2 
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The basic flight condition used was a landing approach along a 3° glide 
slope with flaps set at 0.43633 rad. The approach speed was 77.12 m-sec"^, the 
speed recommended by the manufacturer. 
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PHUGOID STABILITY IN WIND SHEAR 

The phugoid mode may be approximated by assuming that contributions of the 
inertial terms and aerodynamic damping terms to the total applied moment are 
much smaller than those due to changes in u and a. Under these assumptions 
and for = 0, the stability determinant obtained from equation (Al8) reduces 
to 


s - 1 
2 

sin 2 Fq - Xu 

-Xa 

g(cos Fq - bu cos 2 Fq) 



-Zu 

- ^ sin2 Fo 

Uq 

-Za 

-Uq s + g(sin Fq - bu sin 2 Fq) 

= 0 

(Cl) 


-Mu 


0 




This determinant expands to an equation of the form 

a2 OT -L. /i\„ 2 — 


where 


and 


w 




3 + 2^ pcjps + (jL)p — 0 


2?p4)p = -Xu + Mu ^ ^Q (1 - bu cos Vq) 


^Xu - Mu ^j(sin To - bu sin 2ro) - ^Zu - Mu ^j(cos Tq - 


(C2) 


(C3) 


bu cos 2ro) 
(C4) 


As Mu for most aircraft is extremely small oh zero, the expressions for 
2CpWp and tOp^ may be further simplified by neglecting the terms that contain 
Mu* When the Mu terms are omitted, the expressions for 2?pWp and (Op^ 
become 

2Cj^p = -Xu - sib Tod - bu cos Tq) (C5) 



sin Fq - bu 


sin 2Fo) - Zu(cos Fq - bu cos 2Fo)]| 


(C6) 
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Since the characteristic equation (C2) is a quadratic, both 2Cp(»)p and 
0 )p 2 must be positive to guarantee the asymptotic stability of the phugoid; 
otherwise, at least one root of the equation will have a nonnegative real part 
so that the phugoid will be unstable. The curves for 2CpWp = 0 and Wp2 = 0 
as functions of Fg and Oy are shown in figures 13 and 14. These curves 
were calculated by using equations (C5) and (C6) with appropriate data from 
appendix B. This information is combined to determine the regions where equa- 
tion (C2) would be stable or unstable. The regions of stability and instability 
are shown in figure 15 as functions Fg and Oy. Most flight-path angles lie 
in the interval -(tt/18) ^ Fg ^ tt/18, which is the region between the vertical 
lines in figure 15. This figure clearly shows the effect of wind shear dy on 
the phugoid stability in landing (Fg < 0) and level flights (Fg =0). In this 
case, the airplane is stable for all wind shears that have Oy < 1.0. However, 
for climbing flights (Fg > 0), the region of stability is much more restricted 
for negative wind shears. 

That the unstable condition for dy > 0 in figure 15 is solely a function 
of dy and not a function of the airplane parameters can be seen by rewriting 
equation (C6) as 


Uo 


(sin Fg - dy sin 2Fg) - (cos Fg - dy cos 2Fg)_H 


'u 


(C7) 


Since Zy/Xy is the lift-drag ratio of the airplane and is positive with normal 
values between 5 and 16, it can be shown that the sign of 0 )p 2 depends on the 
sign of cos Fg - dy cos^ Fg. To a good degree of accuracy, this means for 
normal flight-path angles that a value of dy > 1.0 will cause the airplane to 
have an unstable phugoid mode. This result is immediate for a small angle 
approximation on Fg. Thus, it is only necessary to determine those combina- 
tions of v^ 1 and Ug that give dy = 1 .0 to determine whether the airplane 
is stable or unstable for positive shear. A curve for estimating this stability 
is shown in figure 4(a). This is a serious type of instability as it arises 
from an interaction of the airplane and its environment and no aerodynamic 
changes to the airplane will correct it. 

The instability for negative shear (dy < 0) (see fig. 15) is not a simple 
function of dy as in the case of positive shear. The effect of negative shear 
on stability can be computed by using equation (C5); however, it is a function 
of the airplane parameters as Xy appears in equation (C5). Generally speak- 
ing, the instability for dy < 0 is a divergent oscillation instead of aperi- 
odic type motion encountered for positive shear (dy > 0). 
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Figure 2.- Updrafts and downdrafts that co-existed with the vertical shear 

shown in figure 1 . 
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Figure 6.- Control-fixed motion of an airplane in wind shear. 
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Figure 7.- The effect of deleting wind acceleration terms in the kinematics 
of the airplane equations when the wind shear is positive, dx = ^u> 

= 0. Motions are control fixed. 



D, distance along approach path, m 

(b) Pitch angle. 


Figure 7.- Continued. 
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Figure 8.- The effect of deleting wind acceleration terms in the kinematics of 
the airplane equations of motion when the wind shear is negative, = Ou, 
Oy = 0. Motions are control fixed. 
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Figure 8.- Concluded. 
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Figure 9.- Block diagram of speed control system. 
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Figure 10.- Comparison of airplane motions in positive wind shear. 
O'p = = 2.0, with and without automatic flight-path and speed 

controls . 
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Figure 10.- Continued. 
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Figure 11.- Comparison of airplane motii 
Ot = Ou “ -2.0, with and without aut< 
speed controls, = 0 . 
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Earth-fixed axes, X^,X2,Xj 


Stability axes, and principal body axes, Xj,X2,X2 

Intermediate axes are not labeled 
and are shown as dashed lines 


Moving axes, X^', X^', xy and principal body axes, x^,X2,x^ 


Figure 12.- Coordinate systems and Euler angles. The order of rotation for 
the Euler angles is 6, and cp. The moving axes translate with air- 

plane and remain parallel to the earth-fixed axes. Positive directions 
are shown. 
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Figure 13.- Combinations of Fq and Ou that make the term 

of equation (C5) zero. 
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Figure 15.- Regions of stability for equation (C2) obtained by combining 

figures 13 and 14. 
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